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We develop a new approach to local nonlinear effects in cosmic microwave background 
anisotropies, and discuss the qualitative features of these effects. New couplings of the baryonic 
velocity to radiation multipoles are found, arising from nonlinear Thomson scattering effects. We 
also find a new nonlinear shear effect on small angular scales. The full set of evolution and constraint 
equations is derived, including the nonlinear generalizations of the radiation multipole hierarchy, and 
of the dynamics of multi-fluids. These equations govern radiation anisotropies in any inhomogeneous 
spacetime, but their main application is to second-order effects in a universe that is close to the 
Friedmann models. Qualitative analysis is given here, and quantitative calculations are taken up in 
further papers. 
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I. INTRODUCTION 

Recent and upcoming advances in observations of the cosmic microwave background (CMB) radiation are fuelling 
the construction of increasingly sophisticated and detailed models to predict the anisotropy on small angular scales. 
Such models require highly specific input in order to produce numerical results, and they involve intricate problems 
of computation. As a complement to such specific predictive models, it is also useful to pursue a more qualitative and 
analytical investigation of CMB anisotropies. A general qualitative analysis does not rely on detailed assumptions 
about the origin of primordial fluctuations, the density parameters of the background, reionization and structure 
formation history, etc. Instead, the aim is to better understand the underlying physical and geometric factors in 
the dynamics of radiation anisotropies, and hopefully to uncover new results and insights. In this paper, we follow 
such an approach, and develop a new analysis of local nonlinear effects in CMB anisotropies. We are able to give a 
physically transparent qualitative analysis of how inhomogeneities and relative motions produce nonlinear effects in 
CMB anisotropies. We derive the nonlinear generalization of Thomson scattering, and we find a new nonlinear shear 
effect on small scales. 

We use a I-I-3 covariant approach (i.e., a “covariant Lagrangian” approach) to CMB anisotropies, based on choice of 
a physically determined 4-velocity vector field u“. This allows us to derive the exact nonlinear equations for physical 
quantities as measured by observers moving with that 4-velocity. Then the nonlinear equations provide a covariant 
basis for investigating second-order effects, as well as for linearizing about a Friedmann-Lemaitre-Robertson-Walker 
(FLRW) background. The basic theoretical ingredients are: (a) the covariant Lagrangian dynamics of Ehlers and 
Ellis and the perturbation theory of Hawking and Ellis and Bruni |Q] which is derived from it; (b) the 1-1-3 
covariant kinetic theory formalism of Ellis, Treciokas and Matravers 1^,^ (which builds on work by Ehlers, Geren and 
Sachs 0, Treciokas and Ellis an d Thorne (H); and (c) the 1-1-3 covariant analysis of temperature anisotropies due 
to Maartens, Ellis and Stoegei^p^. 

The well-developed study of CMB anisotropies is based on the pioneering results in CMB physics (Sachs and Wolfe 
[ pl| , Rees and Sciama [Q, Peebles and Yu [^, Sunyaev and Zeldovich [p^, Grishchuk and Zeldovich and 
others), and on the development of gauge-invariant perturbation theory, particularly by Bardeen and Kodama 
and Sasaki ITtI] (b uilding of the work of Lifshitz ||^). There are comprehensive and detailed models - see e.g. Hu and 
Sugiyama [[19[-j21||, Ma and Bertschinger [^, Seljak et al. |p3|-p^, Durrer and Kahniashvili |^. These provide the 
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basis for sophisticated predictions and comparisons with the observations of recent, current and future satellite and 
ground-based experiments. The hope is that this inter-play between theory and observation (including the large-scale 
galactic distribution and other observations), in the context of inflationary cosmology, will produce accurate values 
for the various parameters that characterize the standard models, thus allowing theorists to discriminate between 
competing models (see for example 

While these papers have provided a near-exhaustive treatment of many of the issues involved in CMB physics, there 
are a number of reasons for pursuing a complementary 1-1-3 covariant approach, as developed in |W 

Firstly, the covariant approach by its very nature incorporates nonlinear effects. This approach starts from the 
inhomogeneous and anisotropic universe, without a priori restrictions on the degree of inhomogeneity and anisotropy, 
and then applies the linearization limit when required. The 1-1-3 covariant equations governing CMB anisotropies are 
thus applicable in fully nonlinear generality. These equations can then be specialized in various ways in addition to a 
standard FLRW-linearization. Second-order effects in an almost-FLRW universe probably form the most important 
possibility, given the increasing accuracy and refinement of observations. The study of CMB anisotropies in homo¬ 
geneous Bianchi universes with large anisotropy is another possibility that flows directly from the general nonlinear 
equations. Such applications will be the subject of future papers in the programme. The current paper is concerned 
with setting up the general dynamical equations and identifying the qualitative nature of nonlinear effects. (The 
general algebraic equations are derived in @|.) 

Secondly, the 1-1-3 covariant approach is based entirely on quantities with a direct and transparent physical and 
geometric interpretation, and the fundamental quantities describing anisotropy and inhomogeneity are all automat¬ 
ically gauge-invariant when a suitable covariant choice of fundamental 4-velocity has been made. As a consequence 
the approach leads to results with unambiguous physical meaning (provided the fundamental 4-velocity field is chosen 
in a physically unique and appropriate way; we discuss the various options below). 

This approach has been developed in the context of density perturbations p|j39|-|^ and gravitational wave per¬ 
turbations PJ55|-|^. (See also for a recent review.) In relation to CMB anisotropies, the covariant Lagrangian 
approach was initiatednby Stoeger, Maartens and Ellis p^ ], who proved the following result: if all comoving observers 
in an expanding universe region measure the anisotropy of the CMB after last scattering to be small, then the universe 
is almost FLRW in that region^ No a priori assumptions are made on the spacetime geometry, or on the source 
and nature of CMB anisotropies, so that this result provides a general theoretical underpinning for CMB analysis 
in perturbed FLRW universes. It effectively constitutes a proof of the stability of the corresponding exact-isotropy 
result of Ehlers, Geren and Sachs [Q. The weak Copernican principle implicit in the assumption that all fundamental 
observers see small anisotropy is in principle partially testable via the Sunyaev-Zeldovich effect (see and references 
therein). The qualitative result was extended into a quantitative set of limits on the anisotropy and inhomogeneity 
of the universe imposed by the observed degree of CMB anisotropy, independently of any assumptions on cosmic 
dynamics or perturbations before recombination |^J3^-^^. 

More recently, this approach to CMB anisotropies in an almost FLRW universe has been extended by Dunsby [Q , 
who derived a 1-1-3 covariant version of the Sachs-Wolfe formula, and by Challinor and Lasenby |35 who performed 
a comprehensive 1-1-3 covariant analysis of the imprint of scalar perturbations on the CMB, confirming the results of 
other approaches from this viewpoint and bringing new insights and clarifications via the covariant approach. In [b 3, 
they also discuss qualitatively the imprint of tensor perturbations on the CMB, in the covariant approach (see 1^ ] 
for quantitative results). 

This paper is closely related to, and partly dependent upon, all of these previous 1-1-3 covariant analyses. It 
extends and generalizes aspects of these papers, using and developing the covariant nonlinear Einstein-Boltzmann- 
hydrodynamic formalism. We analyze the nonlinear dynamics of radiation anisotropies, with the main application 
being second-order effects in an almost FLRW universe. We identify and describe the qualitative features of such 
effects. This lays the basis for a generalization of results on well known second-order effects such as the Rees- 
Sciama and Vishniac effects (see e.g. ||^), and on recent second-order corrections of the Sachs-Wolfe effect [^,Q. 
Developing a quantitative analysis on the basis of the equations and qualitative analysis given here is the subject 
of further work. Ultimately this involves the solution of partial differential equations, which requires in particular a 
choice of coordinates, breaking covariance. However, the 1-1-3 covariant approach means that all the equations and 
variables have a direct and transparent physical meaning. 


^ A 1-1-3 covariant approach to CMB anisotropy was independently outlined by Bonanno and Romano in general terms, 
using a flux-limited diffusion theory, but the detailed implications of small CMB anisotropy were not pursued. 

^ Note the importance of expansion: a static isotropic cosmology with arbitrarily large inhomogeneity can be constructed in 
which all observers see isotropic CMB 
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In Section II, the covariant Lagrangian formalism for relativistic cosmology is briefly summarized. Section III 
develops an exact 1+3 covariant treatment of multi-fluids and their relative velocities, building on In Section 
IV, the covariant Lagrangian approach to kinetic theory is outlined. Section V develops a nonlinear treatment of 
Thomson scattering, which identifies new couplings of the baryonic relative velocity to the radiation multipoles. We 
derive the hierarchy of exact covariant multipole equations which arise from the Boltzmann equation. This section 
uses and generalizes a combination of the results of Ellis et al. on the multipoles of the Boltzmann equation in 
general, Maartens et al. 0 on a covariant description of temperature fluctuations, and Challinor and Lasenby ]^ ] 
on Thomson scattering. The equations constitute a covariant and nonlinear generalization of previous linearized 
treatments. In Section VI, we consider qualitative implications of the nonlinear equations. We identify the role of the 
kinematic quantities in the nonlinear terms, and comment on the implications for second-order effects, which include 
a new nonlinear shear correction to CMB anisotropies on small angular scales. We also give the multipole equations 
for the case where the radiation anisotropy is small, but spacetime anisotropy and inhomogeneity are unrestricted.^ 
Finally, we give the linearized form of the multipole equations, regaining the equations of Challinor and Lasenby 
This provides a covariant Lagrangian version of the more usual metric-based formalism of gauge-invariant 
perturbations (see e.g. In a further paper [^, the linearized equations derived here are expanded in 

covariant scalar modes, and this is used to determine analytic properties of CMB linear anisotropy formation. 

We follow the notation and conventions of , with the improvements and developments introduced by . 

In particular: the units are such that c, SttG and are equal 1; the signature is (—h ++); spacetime indices are 

a, 5, • • • = 0,1, 2, 3; the curvature tensors are W^hcd = —ddy^bc H-? Rab = R'^acb and R = R°‘a, and the Ricci identity 

is V[aVb]itc = \Rabcdu'^] Al denotes the index string 0102 • ■ - a^, and denotes the tensor product • • - 6“^; 

(square) round brackets enclosing indices denote the (anti-)symmetric part, and angled brackets denote the projected 
symmetric and tracefree (PSTF) part (defined below). The spatially projected part of the covariant derivative is 
denoted by D^, following The approximate equality symbol, as in J Ri 0, indicates equality up to first (linear) 

order in an almost-FLRW spacetime. 


II. COVARIANT LAGRANGIAN FORMALISM IN RELATIVISTIG COSMOLOGY 

The Ehlers-Ellis 1+3 formalism |li|,p|J^ is a covariant Lagrangian approach, i.e. every quantity has a natural 
interpretation in terms of observers comoving with the fundamental 4-velocity (where u°‘Ua = — 1). Provided this 
is defined uniquely in an invariant manner, all related quantities have a direct physical or geometric meaning, and 
may in principle be measured in the instantaneous rest space of the comoving fundamental observers. Any coordinate 
system or tetrad can be used when specific calculations are made. These features are a crucial part of the strengths of 
the formalism and of the perturbation theory that is derived from it. We will follow the streamlining and development 
of the formalism given by Maartens the essence of which is to make explicit use of irreducible quantities and 
derivatives, and to develop the identities which these quantities and derivatives obey (see also [|l8y57lJ^ , |^ , 0 ). 

The basic algebraic tensors are: (a) the projector hab = gab + UaUb, where gab is the spacetime metric, which projects 
into the instantaneous rest space of comoving observers; and (b) the projected alternating tensor Sabc = Vabcdu‘^, where 
Vabcd = —\/\g\^^[a^^b^'^c^^d] is the spacetime alternating tensor. Thus 

gabcd ‘^^[a^b]cd i ^abc^ ^ ^b • 

The projected symmetric tracefree (PSTF) parts of vectors and rank-2 tensors are 

V(a) = ha'^Vb, S^ab) = “ \h^''Kb] S^d , 

with higher rank formulas given in 1^ . The skew part of a projected rank-2 tensor is spatially dual to the projected 
vector Sa = ^£abcS^’^'^\ and then any projected rank-2 tensor has the irreducible covariant decomposition 

Rab — "^Rhab + £abcR + R(ab) j 

where R = Rcdh^'^ is the spatial trace. In the 1+3 covariant formalism, all quantities are either scalars, projected vectors 
or PSTF tensors. The equations governing these quantities involve a covariant vector product and its generalization 
to PSTF rank-2 tensors: 


® This case will apply before decoupling, in order to be consistent with the almost-FLRW result quoted above. 
In P J4^ , ^ , ^5|J3^ h is denoted while in ||^,^ it is Va. 
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W^]a = SatcV^W^ , [S, Q]a = SatcS^Q^^^ . 

The covariant derivative Vo defines 1+3 covariant time and spatial derivatives 


= h/h\ ■ ■ ■ h/VdJ" -.../ . 

Note that Dchob = 0 = DijEatc, while hat = 2u(^a'i‘'b) and iabc = 3u[a£bc]dW^- The projected derivative Do further splits 
irreducibly into a 1+3 covariant spatial divergence and curl 

div y = D^Vo , (div S)a = D'’^o6 , 

curl Vo = £obcDV“ , curlS'oh = £cd(aD‘'S'6)‘', 

and a 1+3 covariant spatial distortion 

^{aVb) = 0 ( 0 ^ 6 ) - 3 (divV) hab , 

Note that div curl is not in general zero, for vectors or rank-2 tensors (see [l^ , ^8|j^ , |5^ for the relevant formulas). 
The covariant irreducible decompositions of the derivatives of scalars, vectors and rank-2 tensors are given in exact 
(nonlinear) form by 

VoV' = -'0Mo + DaV', (1) 

VbVo = -Ub {V(o> + A,V-Ua] + Uo {^e^b + (TbcV^ + [w, V]b} 

+ i (divV) hab - 5£obcCuriy‘= + D(oV),) , (2) 

^cSab — ■^‘^(ob) T ^‘^{a^b)dA + 2u^o {3®‘^b)c T ^b) (^^cd ^cde^ )} 

+ I {divS)^^hi,)a - ^£dc(aCnr\Sb)‘^ + J^{aSbc) ■ (3) 

The algebraic correction terms in equations i) and (^ arise from the relative motion of comoving observers, as 
encoded in the kinematic quantities: the expansion 0 = D“ua, the 4-acceleration Aa = Ua = A(^a)^ the vorticity^ 
uja = —^curlua, and the shear dab = D^a^^b)- Thus, by Eq. (||) 

^b^a — AaUb + ^Qhab + ^abc^ + ^ab ■ 

The irreducible parts of the Ricci identities produce commutation identities for the irreducible derivative operators. 
In the simplest case of scalars: 


CUrlDaV^ = £QbcD['’D^V = -S'f/’Wa , (4) 

DaV' - ha^ (DbV’) = -V'^a + |0DoV' + + [bJ, ^tp]a ■ (5) 

Identity (^ reflects the crucial relation of vorticity to non-integrability; non-zero tOa implies that there are no constant¬ 
time 3-surfaces everywhere orthogonal to since the instantaneous rest spaces cannot be patched together smoothly.^ 
Identity (|) is the key to deriving evolution equations for spatial gradients, which covariantly characterize inhomo¬ 
geneity 1^ I . Further identities are given in @@,§,1311. 

The kinematic quantities govern the relative motion of neighboring fundamental world-lines, and describe the 
universal expansion and its local anisotropies. The dynamic quantities describe the sources of the gravitational field, 
and directly determine the Ricci curvature locally via Einstein’s field equations. They are the (total) energy density 
p = TabU°'u^, isotropic pressure p = \habT°‘^, energy flux qa = —Ti^a)b'u^^ and anisotropic stress TTab = T(ab >5 where Tab 
is the total energy-momentum tensor. The locally free gravitational field, i.e. the part of the spacetime curvature not 


® The vorticity tensor ujab = Sabcb)‘^ is often used, but we prefer to use the irreducible vector uJa- The sign conventions, 
following l^j^, are such that in the Newtonian limit, a; = — x a. Note that D*’a;a6 = curlwa. 

® In this case, which has no Newtonian counterpart, the Da operator is not intrinsic to a 3-surface, but it is still a well-defined 
spatial projection of Va in each instantaneous rest space. 
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directly determined locally by dynamic sources, is given by the Weyl tensor Cabcd- This splits irreducibly into the 
gravito-electric and gravito-magnetic fields 

^ah — Cachd'd’ U — ^[ah) 5 ^ab — 2^acdC he'^ — ^{ab) j 

which provide a covariant Lagrangian description of tidal forces and gravitational radiation. 

An FLRW (background) universe, with its unique preferred 4-velocity is covariantly characterized as follows: 

dynamics: = 0 = DaP, Qa = 0, T^ab = 0; 

kinematics: Da0 = 0, Aq = 0 = uja, (Jab = 0; 
gravito-electric/magnetic field: Eab = 0 = Hab- 

The Hubble rate is i/ = ^0 = d/a, where a{t) is the scale factor and t is cosmic proper time. In spatially 
homogeneous but anisotropic universes (Bianchi and Kantowski-Sachs models), the quantities qa, T^ab, o'ab, Eab and 
Elab in the preceding list may be non-zero. 

The Ricci identity for and the Bianchi identities V^Cabcd = V[a(—i?Wc + g^5h]c) produce the fundamental evo¬ 
lution and constraint equations governing the above covariant quantities Hi- Einstein’s equations are incorporated]] 
via the algebraic replacement of the Ricci tensor Rab by Tab — ^Tc'^gat- These equations, in exact (nonlinear) form 
and for a general source of the gravitational field, are : 

Evolution: 


P+{P + P)0 + div g = -2A“ga - cr“*'7rab , (6) 

0 -I- i0^ -I- i(p-(- 3p) - divA = -(Tab(^°‘'' + 2u;aw“ -I- A^A" , (7) 

9(a) + |09a + (p + p)Aa + T>aP + (div 7r)a = -(Jabq^ +[uj,q]a- A^TTab , (8) 

W(a) + §0Wa + ^CUrl Aa = CTafew'’ , (9) 

^{ah) T 30 f^afe “f Eab ' 2 '^ah B^qA^^ ^c{a^b) ^{a^b) T A(^aAb) , (^0) 

E{ab) + QEab - CnvlHab + \{p+p)<yab + {ab) + \P(^adb) + l^T^ab 

= ~ A(^a<lb) + ‘^A'^ecd{aHb)‘^ + 3cra{aEb)‘^ 

^ ^cd{aEb) 2 *^ {a'^b)c 2 ^ ^cd{a'^b) 7 (H) 

E{ab) T ^Eab T CUrl TTab c{aEb) ^ ^cd{aEb) 

- ‘^'A'^ScdiaEb)'^ - |^^(a9fc) + ^Cr‘'{a£b)cdq‘^ ■ (12) 

Constraint: 

dlYUJ = A°'UJa , (13) 

(div(T)a - curlwa - |Da0 + qa = -2[u;, A]a , (14) 

curl (Tab + D^ai^b) - Hab = -2A(aWb) , (15) 

(divE;)Q -I- i(div7r)Q - |DaP-f ^Qqa = [o',H]a - iHabu/’ + \aabq^ - |[w,9]a , (16) 

(div77)a + ^curlga - (p + p)Wa = -[u,E]a - ^[<J,TT]a + SE’abw'' - ^TTabUj'’ . (17) 


If the universe is close to an FLRW model, then quantities that vanish in the FLRW limit are 0(e), where e is 
a dimensionless smallness parameter, and the quantities are suitably normalized (e.g. \/o//ba^/H < e, etc.). The 
above equations are covariantly and gauge-invariantly linearized |Q| by dropping all terms O(e^), and by replacing 
scalar coefficients of 0(e) terms by their background values. This linearization reduces all the right hand sides of the 
evolution and constraint equations to zero. 


III. 1+3 COVARIANT NONLINEAR ANALYSIS OF MULTI-FLUIDS 

The formalism described above applies for any covariant choice of If the physics picks out only one then 
that becomes the natural and obvious 4-velocity to use. In a complex multi-fluid situation, however, there are various 


^ Note that one constraint Einstein equation is not explicitly contained in this set - see 
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possible choices. The different particle species in cosmology will each have distinct 4-velocities; we could choose 
any of these as the fundamental frame, and other choices such as the centre of mass frame are also possible. This 
allows a variety of covariant choices of 4-velocities, each leading to a slightly different 1-1-3 covariant description. One 
can regard a choice between these different possibilities as a partial gauge-fixing (but determined in a covariant and 
physical way). Any differences between such 4-velocities will be 0{e) in the almost-FLRW case and will disappear 
in the FLRW limit,^ as is required in a consistent 1-1-3 covariant and gauge-invariant linearization about an FLRW 
model (see I® for further discussion). 

In addition to the issue of linearization, one can also ask more generally what the impact of a change of fundamental 
frame is on the kinematic, dynamic and gravito-electric/magnetic quantities. If an initial choice is replaced by a 
new choice then 


= 7 (u“-I-n“) where VaU°'= 0 , j = {I — VaV^-) , (18) 

where is the (covariant) velocity of the new frame relative to the original frame. The exact transformations of all 
relevant quantities are given in the appendix, and are taken from [ ^ . To linear order, the transformations take the 
form: 


0«0-|-divV, Aa ^ Aa + Va + HVa , 

UJa ^ UJa - ^CuAVa , dah « (Tab -f D, 

p, P^P, qa^ qa- ip + p)Va , ^ab « T^ab , 
^ab ~ ^ab ^ ^ab ~ ■ 


Suppose now that a choice of fundamental frame has been made. (For the purposes of this paper, we will not 
need to specify such a choice.) Then we need to consider the velocities of each species which source the gravitational 
field, relative to the fundamental frame. If the 4-velocities are close, i.e. if the frames are in non-relativistic relative 
motion, then 0{v^) terms may be dropped from the equations, except if we include nonlinear kinematic, dynamic and 
gravito-electric/magnetic effects, in which case, for consistency, we must retain 0(e°u^) terms such as pv^, which are 
of the same order of magnitude in general as O(e^) terms. (See ^^.) If the universe is close to FLRW, then 0{e°v'^) 
terms may be neglected, together with 0{ev) and O(e^) terms. 

In summary, there are two different linearizations: 

(a) linearizing in relative velocities (i.e. assuming all species have nonrelativistic bulk motion relative to the 
fundamental frame), without linearizing in the kinematic, dynamic and gravito-electric/magnetic quantities that 
covariantly characterize the spacetime; 

(b) FLRW-linearization, which implies the special case of (a) obtained by also linearizing in the kinematic, dynamic 
and gravito-electric/magnetic quantities. 

Clearly (a) is more general, and we can take it to be the physically relevant nonlinear regime, i.e. the case where 
only nonrelativistic average velocities^ are considered, but no other assumptions are made on the physical or geometric 
quantities. In case (a), no restrictions are imposed on non-velocity terms, and we neglect only terms 0{ev'^,v^). In 
case (b), we neglect terms 0{e‘^, ev,v‘^). Covariant second-order effects against an FLRW background are included 
within (a), when we neglect terms O(e^). (Note that gauge-invariance is a far more subtle problem at second order 
than at first order: see Bruni et al. [7^].) 

The dynamic quantities in the evolution and constraint equations (i)-® are the total quantities, with contribu¬ 
tions from all dynamically significant particle species. Thus 


T ab \ '' rriab a b . ^ucib , rt (a b) , _ab 

= y. = pu u pfi + 2q^ u ' TT , 
1 

t;^ = + p,hf + 2q‘^'^y + , 


(19) 

( 20 ) 


where / labels the species. We include radiation photons (/ = R), baryonic matter (/ = B) modelled as a perfect 
fluid, cold dark matter (/ = C) modelled as dust over the era of interest for CMB anisotropies, neutrinos (/ = N) 


® A similar situation occurs in relativistic thermodynamics, where suitable 4-velocities are close to the equilibrium 4-velocity, 
and hence to each other | [n| . 

® Of course, this implies no restrictions on the velocities of individual particles within any species. 
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(assumed to be massless), and a cosmological constant (/ = 1^).^ Note that the dynamic quantities p,, • • ■ in equation 
(^) are as measured in the /-frame, whose 4-velocity is given by 

= 7/ (w“ + 1^“) : 'V^Ua = 0 . (21) 

Thus we have 


^ah 


Pa = 0 = C = 

Pr sPr ’’ Pn 3 Av > 




ab 


( 22 ) 

(23) 


where we have chosen the unique 4-velocity in the cold dark matter and baryonic cases which follows from modelling 
these fluids as perfect. The cosmological constant is characterized by 

Pv = -A. = -A, = 0 = 7^“^ = 0. 


The conservation equations for the species are best given in the overall it“-frame, in terms of the velocities of 
species / relative to this frame. Furthermore, the evolution and constraint equations of Section II are all given in 
terms of the M“-frame. Thus we need the expressions for the partial dynamic quantities as measured in the overall 
frame. The velocity formula inverse to equation (|^) is 


= 7i « + <“) , = -7/ « + 


(24) 


= 0, and = v^Vja- Using this relation together with the general transformation equations 

(A8)-( All ), or directly from the above equations, we find the following exact (nonlinear) equations for the dynamic 
quantities of species / as measured in the overall u“-frame: 


Pj = Pi + {P, +Pi) + 27/gfUa + : 

p* =Pi + l {7,V (a +Pr) + > 

+ {p, +pj< 

+ {(T - 1 ) 9 “ - 7/A^Ab’^“ + if'vf (Pi +PrX + , 

^*ab ^ ^ab ^ 

+ {-iKScVjdh^^ + f {p, + + 27,n/“g;’>} . 


(25) 

(26) 

(27) 

(28) 


These expressions are the nonlinear generalization of well-known linearized results (see e.g. p5||fl| ]). FLRW lineariza¬ 
tion implies that <C 1 for each /, and we neglect all terms which are 0{v^) or 0{eVj). This removes all terms in 
braces, dramatically simplifying the expressions: 


P* ~ Pi , P* ~ Pi , 9 /“ ~ 9 “ + (a + a)a“ ’ ~ 


^ah 


To linear order, there is no difference in the dynamic quantities when measured in the /-frame or the fundamental 
frame, apart from a simple velocity correction to the energy flux. But in the general nonlinear case, this is no longer 
true. 

The total dynamic quantities are simply given by 


III I 

Note that the equations (^^-(|2|) have been written to make clear the linear parts, so that the irreducible nature is 
not explicit. Irreducibility (in the M“-frame) is revealed on using the relations 


A more general treatment, incorporating all the sources which are currently believed to be potentially significant, would 
also include a dynamic scalar field that survives after inflation (“quintessence”), and hot dark matter in the form of massive 
neutrinos (see for a survey with further references). Our main aim is not a detailed and comprehensive model with numerical 
predictions, but a qualitative discussion focusing on the underlying dynamic and geometric effects at nonlinear and linear level 
that are brought out clearly by a 1-1-3 covariant approach. In principle our approach is readily generalized to include other 
sources of the gravitational field. 
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= h\q^ = - q^bU^ , 

T^ia)ib} ^ ^ ^ab _ + n^\,V,dU^u’^ . 

The exact equations show in detail the specihc couplings and contributions of all partial dynamic quantities in the 
total quantities. For example, it is clear that in spatially homogeneous but anisotropic models, the partial energy 
fluxes contribute to the total energy density, pressure and anisotropic stress at first order in the velocities Vj , while 
the partial anisotropic stresses 7 r“^ contribute to the total energy flux at first order in Vj. 

The total and partial 4-velocities define corresponding number 4-currents: 


= nu“ +f = Y^ iV/ , 

I 




(29) 


where n and Uj are the number densities, and are the number fluxes, and jaU°' = 0 = j,aU^- It follows that 

n = ^ n; = ^ n, -h ^ {( 7 , - l)n, -F j^^v,a} , (30) 

III 

f = Er = E +^^ 0 +E {( 0 ^ -’ (31) 


where the starred quantities are as measured in the it“-frame. Linearization removes the terms in braces, regaining 
the expressions in [ ph|j7l|] . 

Four-velocities may be chosen in a number of covariant and physical ways. The main choices are |75 0 : (a) the 
energy (Landau-Lifshitz) frame, defined by vanishing energy flux, and (b) the particle (Eckart) frame, defined by 
vanishing particle number flux. For a given single fluid, these frames coincide in equilibrium, but in general they are 
different. For each partial an y chan ge in choice leads to transformations in the partial dynamic quantities, 

that are given by equations (A 8 )-( All ) in the appendix. For the fundamen tal a change in choice leads in addition 
to transformations of the kine m atic quantities, given by equations (A4)-(A7), and of the gravito-electric/magnetic 
field, given by equations ( |A12| )-( [A1^ ). 

A convenient choice for each partial four-velocity is the energy frame, i.e. q^ = 0 for each / (this is the obvious 
choice in the cases I = C, B). As measured in the fundamental frame, the partial energy fluxes do not vanish, i.e. 

7 ^ 0 , and the total energy flux is given by 


9 “ = E [(ft + ft) + 0{€V^,V^)] . 


(32) 


With this choice, using the above equations, we find the following expressions for the dynamic quantities of matter 
as measured in the fundamental frame. For cold dark matter: 


For baryonic matter: 


Pc = ^cPc > Pc 
C = IcPc^c ’ 


= - 


pI = 7b (1 + Ps, pI = [wb + + W^)] Pb , 

C = 7b(1 + Wj,)psV^ , 7r;“*' = 7^(1 -h Wb)Pbi4“^'b^ > 


(33) 

(34) 


(35) 

(36) 


where = Vb!Pb- I'^ Hi® case of radiation and neutrinos, we will evaluate the dynamic quantities relative to the 

M“-frame directly via kinetic theory, in the next section. 

The total energy-momentum tensor is conserved, i.e. = 0, which is equivalent to the evolution equations (|^) 

and (^). The partial energy-momentum tensors obey 

VfcT/'’ = J/ = C/;u“ -f , (37) 


where B* is the rate of energy density transfer to species I as measured in the it“-frame, and is the 

rate of momentum density transfer to species /, as measured in the u“-frame. Cold dark matter and neutrinos are 
decoupled during the period of relevance for CMB anisotropies, while radiation and baryons are coupled through 
Thomson scattering. Thus 
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where the Thomson rates are 


-/“ = 0 = j;, 4“ = -j; = C4w“ + m; , 


(38) 


Uj, - TL^ (Jrj. 

= ^sCTt 


[Ip>b 


-<C^Ba) +0{evl,vl), 

-C + CV)+0(6.,^ 



(39) 

(40) 


as given by Eq. (^), derived in Section V. Here is the free electron number density, and oy is the Thomson 
cross-section. Note that to linear order, there is no energy transfer, i.e. C4 ~ 0. 

Using equations ([3^)-(|^) in (^, we find that for cold dark matter 


Pc + + Pcdivib = - (Pc^^c) - 

- v^^aPc - ‘^Pc^aV^ + 0{ev^,v ^), (41) 

+ -^c^bv^ + 0{ev^,v^), (42) 

and for baryonic matter 

Pb + 0(1 + 'Ws)ps + (1 + Ws)pgdivvs 

= - [{l+W^)p^V^] - ^V^e{l+W^)p^ -<Da [(1 +Wb)^,] 

- 2(1 -f w^)p^AaV^ - (|p>b “ + 0(eu^, uj), (43) 

(1 -I- u;b)u“ + (I - c^) + (1 + Wg)A°‘ 

+ Pb + pB^n^<^T (P>B - O 

= {l + w^)Abvlu^ - {l + w^)a^bvl + {1 + Wj,)[uj,v^]‘^ 

- (1 + + clil + w,){diyv,X 

-PB^’^BOT<“^^'B6 + 0(eu^,uj), (44) 


where j L 

those given in IS 


(this equals the adiabatic sound speed only to linear order). These conservation equations generalize 
I to the nonlinear case. FLRW linearization reduces the right hand sides of these equations to zero, 
dramatically simplifying the equations. The conservation equations for the massless species (radiation and neutrinos) 
are given below. Note from Eq. ( ^^ that if the cold dark matter frame is chosen as the fundamental frame^hen the 


4-acceleration vanishes, i.e. = 0 implies Aa = 0. This is the choice of fundamental frame advocated in 361. 


IV. COVARIANT LAGRANGIAN KINETIC THEORY 

Relativistic kinetic theory (see e.g. |7^-|7^) provides a self-consistent microscopically based treatment where there 
is a natural unifying framework in which to deal with a gas of particles in circumstances ranging from hydrodynamic 
to free-streaming behavior. The photon gas undergoes a transition from hydrodynamic tight coupling with matter, 
through the process of decoupling from matter, to non-hydrodynamic free streaming. This transition is characterized 
by the evolution of the photon mean free path from effectively zero to effectively infinity. The range of behavior can 
appropriately be described by kinetic theory with Thomson scattering and the baryonic matter with which 

radiation interacts can reasonably be described hydrodynamically during these times. (The basic physics of radiation 
and matter and density perturbations in cosmology was developed in the works of Sachs and Wolfe ||^, Silk [ pl| , 
Peebles and Yu ||^, Weinberg |^, and others.) 

In the covariant Lagrangian approach of (see also |^,^), the photon 4-momentum (where p“pa = 0) is split as 

= E(m“ -I- e“), e“ea = 1, eAua = 0 , (45) 

where E = —UaP°' is the energy and e“ = jE is the direction, as measured by a comoving (fundamental) observer. 
Then the photon distribution function is decomposed into covariant harmonics via the expansion 

f{x,p) = f{x,E,e) = F + Fae°’ + Fabe‘^e^-\ -= ^ Ea, (x, (46) 

e>o 
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where and provides a representation of the rotation group [^. The covariant multipoles 

are irreducible since they are PSTF, i.e. 

Fa-.-b = F(^a- b) ^ Fa-.-b = F(^a- b) , Fa-■-bU^ = Q = Fa-■ ■bch^‘^ ■ 

They encode the anisotropy structure of the distribution in the same way as the usual spherical harmonic expansion 

i>0 m——£ 


but here (a) the Fa^ are covariant, and thus independent of any choice of coordinates in momentum space, unlike the 
/™; (b) Fai is a rank-£ tensor field on spacetime for each fixed E, and directly determines the £-multipole of radiation 
anisotropy after integration over E. The multipoles can be recovered from the distribution function via 1^,^^ 


Fai = A; 


y/(xjif, e)e(A^)dn, with Ae = An- 


(£!) 22 ^ 
( 2 ^+ 1 )! 

where df2 = d^e is a solid angle in momentum space. A further useful identity is |^] 

0 i odd, 


(47) 


e^‘dn = 


47r 


[ hi^i‘>.2l^a3a4 _ ,f^at-iai) g 

The first 3 multipoles arise from the radiation energy-momentum tensor, which is 
Tf{x) = I pV/(^,P)d"p = p;u‘^u^ + + ■ 


(48) 


.*ab 


where d^p = EdEd^ is the covariant volume element on the future null cone at event x. It follows that the dynamic 
quantities of the radiation (in the M“-frame) are: 


p*=47r/ E^FdE,q*^ = 


Att 


> o_ poo 

E^E^dE, = 




(49) 


From now on, we drop the asterisks from the radiation dynamic quantities relative to the fundamental frame, since 
we do not need to relate them to their values in the radiation frame. 

We extend these dynamic quantities to all multipole orders by defining^ ^ 


Ilai -aj — / E^Ea^-.-agdE , 


(50) 


so that n = p„/47r, n“ = 3g“/47r and n“^ = IStt^VStt. 
The Boltzmann equation is 


dt! dx°‘ ^ dp°- ’ 


(51) 


where p“ = dx°‘/dv and C[f] is the collision term, which determines the rate of change of / due to emission, absorption 
and scattering processes. This term is also decomposed into covariant harmonics: 


E[f] = X! E)e^‘ =b + bae°‘ + babe°’e^ H-, 

l>0 


(52) 


where the multipoles bAf = b(^Ae) encode covariant irreducible properties of the particle interactions. Then the 
Boltzmann equation is equivalent to an infinite hierarchy of covariant multipole equations 


Because photons are massless, we do not need the complexity of the moment definitions used in Q]. In [Q, is used, 
where = AeTlAf From now on, all energy integrals will be understood to be over the range 0 < < oo. 
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Lai{x,E) = bAti^ArnK^^E) 


where La^ are the multipoles of d//di;, and will be given in the next section. These multipole equations are tensor 
field equations on spacetime for each value of the photon energy E (bu t note that energy changes along each photon 
path). Given the solutions Fai{x,E) of the equations, the relation (ff^ ) then determines the full photon distribution 
f{x, E, e) as a scalar field over phase space. 

Over the period of importance for CMB anisotropies, i.e. considerably after electron-positron annihilation, the 
average photon energy is much less than the electron rest mass and the electron thermal energy may be neglected, 
so that the Compton interaction between photons and electrons (the dominant interaction between radiation and 
matter) may reasonably be described in the Thomson limit. (See [Q for refinements.) We will also neglect the effects 
of polarization (see e.g. Q). For Thomson scattering 

C[f] = cr^n^E^ [f{x,p) - f{x,p)] , (53) 


where E^ = —paU^ is the photon energy relative to the baryonic (i.e. baryon-electron) frame u“, and f{x,p) determines 
the number of photons scattered into the phase space volume element at ix,p). The differential Thomson cross-section 
is proportional to 1 -I- cos^ a, where a is the angle between initial and final photon directions in the baryonic frame. 
Thus cos a = where is the initial and e“ is the final direction, so that 


p'“ = F;,«+e-) , = + , 

where we have used which follows since the scattering is elastic. Here is given by Eq. (^), where is 

the velocity of the baryonic frame relative to the fundamental frame with v^Ua = 0. Then / is given by |^J7^ 


fix,p) = 


1071 


dH;. 


(54) 


f{x,p') [l + (eXa)^ 

The exact forms of the photon energy and direction in the baryonic frame follow on using equations and ([^): 

(55) 

(56) 


E, = E^, (1 - <ea) , 

^ [e“ -f yj {v^et - v^) -f yj (v^eb - l) <] . 




7b( 1 -^|ec) 


Anisotropic scattering will source polarization, and small errors are introduced by assuming that the radiation 
remains unpolarized . A fully consistent and general treatment requires the incorporation of polarization. However, 
for simplicity, and in line with many previous treatments, we will neglect polarization effects. 


V. THE NONLINEAR MULTIPOLE HIERARCHY 

The full Boltzmann equation in photon phase space contains more information than necessary to analyze radiation 
anisotropies in an inhomogeneous universe. For that purpose, when the radiation is close to black-body we do 
not require the full spectral behaviour of the distribution multipoles, but only the energy-integrated multipoles. The 
monopole leads to the average temperature, while the higher order multipoles determine the temperature fluctuations. 
The 1-1-3 covariant and gauge-invariant definition of the average temperature T is given by [l^ 

p^{x) j E^F{x,E)AE = rT{xf , (57) 

where r is the radiation constant. If / is close to a Planck distribution, then T is the thermal black-body average 
temperature. But note that no notion of background temperature is involved in this definition. There is an all-sky 
average implied in (^). Fluctuations across the sky are measured by integrating the higher multipoles (a precise 
definition is given below), i.e. the fluctuations are determined by the liai - at > 1) defined in Eq. ([^). 

The form of C[f] shows that covariant equations for the temperature fluctuations arise from decomposing the 
energy-integrated Boltzmann equation 


/ 


E^^dE 

dv 


E^C[f]dE 


(58) 
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into 1+3 covariant multipoles. We begin with the right hand side, which requires the covariant form of the Thomson 
scattering term (]^). Since the baryonic frame will move nonrelativistically relative to the fundamental frame in all 
cases of physical interest, it is sufficient to linearize only in , and not in the other quantities. Thus we drop terms in 
0(ev^,v^) but do not neglect terms that are 0(e°u|, erig) or O(e^) relative to the FLRW limiting background. In other 
words, we make no restrictions on the geometric and physical quantities that covariantly characterize the spacetime, 
apart from assuming a nonrelativistic relative average velocity for matter. The resulting expression will in particular 
be applicable for covariant second-order effects in FLRW backgrounds (recognising that polarization effects should be 
included for a complete treatment), or for first-order effects in Bianchi backgrounds. 

For brevity, we will use the notation 

0[3]^0(eug^uj), 

noting that this does not imply any second-order restriction on the dynamic, kinematic and gravito-electric/magnetic 
quantities. It follows from equations (^) and that 

47r J jE^dE^ = (pg)g + |«'’)gegaegf,, (59) 


where the dynamic radiation quantities are evaluated in the baryonic frame. This approach relies on the frame- 
transformations given in the appendix, and allows us to evaluate the Thomson scattering integral more directly and 
clearly tha n oth er approaches. In the process, we are also generalizing to include nonlinear effects. We use equations 
(AS) and (All) to transform back to the fundamental frameiQ 


(Pb)b = Ai [1 + I^'b] - 2g“r;Ba + C>[3], 

K')b = K' + + 0[3] 


Now 


E^C[f]dE = 1 + - 


2 


[l-<e, + iu2] J fE^dE^ 


E!fdE, 


0[3]. 


(60) 


In addition, we need the following identity, valid for any projected vector 

+ [Fv, + e“ 

+ [F(aVb) + jFabcV^^] + • ■ • 


= E 


^>0 ■- 


F{At-iVai) + 


e + 1 

2 £ + 3 


FA.aV'^ 


Ml) 


(61) 


(Here and subsequent^, we use the convention that = 0 for ^ < 0.) This identity may be proved using Eq. (|4 
and the identity (see H, p. 470): 


V{bSA,) = V(bSA,) - V^S^(A,_^K,b) where Sa, = S^a,) ■ 

Using the above equations, we find that^ 

An J E^C[f]dE = [^PmI - g>fla] 


(62) 


- ngCq, [3g“ - 4pgUg - 37r“^UB6] 






— 127 ^ 77 “^'^?; - 3o 


^{a^h) 


_ l^^nabcd 


Vsd 


^{aSb^c) + • ■ ■ + 0[3] . 


(63) 


As noted in Section III, we retain the 0(Wg) term in (pg)g since pg is zero-order. 
A. Challinor has independently derived the same resnlt ||38|. 
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Now it is clear from equations (^) and that the first four multipoles are affected by Thomson scattering differently 
than the higher multipoles. This is confirmed by the form of equation (|^). Defining the energy-integrated scattering 
multipoles 


Ka, 


EHaAE, 


we find from Eq. 


|) that 

K = [lUv^ - in^,] + 0[3], 

= -n^a^ [n“ - 4 nu“ - - 


■ C^[3], 




T^abc _ _ 

K = —n^CFj. 




■ 0[3], 


jjabc _ 3n{ab^c> _ 


0[3] 


and, for £ > 3: 


(64) 

(65) 

( 66 ) 

(67) 


K ‘ = -UgOj, 




£ + 1 
2£-£3 


n 


-A^a 


+ 0[3] 


( 68 ) 


Equations @-(11) are a nonlinear generalization of the results given by Challinor and Lasenby They 

show the new coupling of baryonic bulk velocity to the radiation multipoles, arising from loeal nonlinear effects in 
Thomson scattering. If we linearize fully, i.e. neglect all terms containing except the term in the dipole 
K°', which is first-order, then our equations reduce to those in [Q. The generalized nonlinear equations apply to 
the analysis of second-order effects on an FLRW background, to first-order effects on a spatially homogeneous but 
anisotropic background, and more generally, to any situation where the baryonic frame is non-relativistic relative to 
the fundamental u“-frame. 

Next we require the multipoles of d//du. These can be read directly from the general expressions first derived in 
iQ, which are exact, 1-1-3 covariant and also include the case of massive particles. For clarity and completeness, we 
outline an alternative, 1-1-3 covariant derivation (the derivation in [|j uses tetrads). We require the identity |@,|| 


du 


— E^ [^0 + Aae^" + CTahC^e^] , 


(69) 


which follows directly from E = —p°-Ua, p^VbP°‘ = 0 and VbUa = —AaUb + D^Ua- Then 


— [Ea,...aAx,E)e° 


= ^[E-^Ea,...a,A,E)p‘^^---p‘^^] 


= E { [i0 + AbC^ + abccA^] {iEa,...a, - EFA..J ■ • 
+ (u“^ + + e“0 [K-a^ + e^ybEa,...a,] } 


where a prime denotes d/dE. The first term is readily put into irreducible PSTF form using the identity (^) with 
Va = Aa, and its extension to the case when 14 is replaced by a rank-2 PSTF tensor Wab (see §, p. 470), with 
Wab = ofab- In the second term, when the round brackets are expanded, only those terms with at most one u°‘'' survive, 
and 


u°-Fa... = -A^Fa..., AV^Fb... = - (i 0 /i“*' -f Eb... . 

Thus the covariant multipoles bA^ of d//du are 


(£ + 1 ) 
( 2 £T^ 

’+ 1 ) z.- 


E ^bAi — Aai) - |0-£'£^Af + ^{aeEAe_i) + 




(2£-f 3) 




at) 


_p,,b^ P e (£-I- 1)(£-I- 2) \ T7 , l' afc 

£ £bc(at At-t) (2£-h3)(2£-f 5)^ ^ 


2 £ 


(2£-f 3) 


E 


- 1/2 


E^/^Fb^AtJ'^at)” - E^-^ [E^-^F^At-.Uat-^at) • 


(70) 
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This regains the result of Q [equation (4.12)] in the massless case, with minor corrections. The form given here 
benefits from the streamlined version of the 1+3 covariant formalism. We reiterate that this result is exact and holds 
for any photon or (massless) neutrino distribution in any spacetime. We now multiply Eq. ( fz^ ) by and integrate 
over all energies, using integration by parts and the fact that E’^Fa-.. —> 0 as i? ^ oo for any positive n. We obtain 
the multipole equations that determine the brightness multipoles Ha^ ■ 


Ka, 


n(Af) + + 


(^+ 1 ) 

+ 3) 




^ ^^bAt + (« + - 

(£-!)(£+!)(£ +2) , 5£ , 

(2^ + 3)(2^ + 5) (2£ + 3) 


i-UJ £bc(a£^A£^i) 




i^At- 2 ) ■ 


(71) 


Once again, this is an exact result, and it holds also for any collision term, i.e. any Kai- For decoupled neutrinos, 
we have = 0 in this equation. For photons undergoing Thomson scattering, the left hand side of Eq. (0) is 
given by Eq. (| 6 ^, which is exact in the kinematic and dynamic quantities, but first order in the relative baryonic 
velocity. The equations (|^) a nd ( 0 ) thus constitute a nonlinear generalization of the FLRW-linearized case given 
by Challinor and Lasenby||^. 

These equations describe evolution along the timelike world-lines of fundamental observers, not along the lightlike 
geodesics of photon motion. The timelike integration is related to light cone integrations by making homogeneity 
assumptions about the distribution of matter in (spacelike) surfaces of constant time, as is discussed in Qj. 

The monopole and dipole of equation 0 ) give the evolution equations of energy and momentum density: 

K = n + len + , (72) 

= n<“> + | 0 n“ + D“n + 

+|^n“'’ + 4 nA“ - [w, n]“ + a^'^Ub. (73) 

In the case of neutrinos, K,^ = 0 = 77“, these express the conservation of energy and momentum]^ 


Pn s^Pn ~ ^ab'^j 


ab 




= -[^, - Abirf . 


FLRW-linearization reduces the right hand sides to zero. For photons, K and 77“ are given by equations 
(1^), and determine the Thomson rates of transfer in equations (39) and 


TL = 47r77, 


47r 

M“ = —77“ 
^ 3 


(74) 

(75) 
and 

(76) 


Finally, we return to the definition of temperature anisotropies. As noted above, these are determined by the 
Ua/:- Generalizing the linearized 1+3 covariant approach in ||^, we define the temperature fluctuation t{x, e) via the 
directional bolometric brightness: 


T{x) [l + r(a:,e)] 


47r 

r 


1/4 


E^f(x,E,e)dE 


(77) 


This is a 1+3 covariant and gauge-invariant definition which is also exact. We can rewrite it explicitly in terms of the 
Aai- 


T{x,e) = 


1 + 



- 1 = Tae“ + Tabe°-e^ -\ -. 


(78) 


As in the photon case, we omit the asterisks on the neutrino dynamic quantities, since we do not require their values in the 
neutrino frame. 
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In principle, we can extract the irreducible PSTF temperature fluctuation multipoles by using the inversion in Eq. 

( 13 ): 


taAx) = ^ 


T(x,e)e^At}dfl. 


(79) 


In the almost-FLRW case, when r is 0(e), we regain from Eq. (|78|) the linearized definition given in |p^: 


TAt 


PrJ 


(80) 


where £ > 1. In particular, the dipole and quadrupole are 


, 157r, 

—— and «- - 

dPn ‘^Pr 


ah 


(81) 


VI. QUALITATIVE IMPLICATIONS OF THE NONLINEAR DYNAMICAL EFFECTS 


In Section II, we gave the nonlinear evolution and constraint equations governing the kinematic, total dynamic and 
gravito-electric/magnetic quantities - see equations i)-( 0 )- In these equations, the total dynamic quantities are, 
using the results of Section III: 


P — Pr + Pm + (l + v^) Pc + \l + (1 + Wg)vl\ pg + K +0[^\, (82) 

P=\Pr + \Pn + I^Ipc + [wb + i(l + )vl] Pfl - A + C>[3], (83) 

+ (1 + 'Wb)PbVb + i (84) 

^ ^ ^ab ^ ^ ^ W,)p,V^^V^'> + 0 [ 3 ] . ( 85 ) 


The conservation equations for matter were given in Section III - see equations (^l|)-(^4|). For neutrinos, the equations 
were given in Section V - see equations (|^) and (!D- For photons, the equations follow from the results of Section 
V as: 


+ |0g“ + + D(,7r 

I b r, . ^ I 


,ab 

ab 


itPR^B - €^Bo) + 0[3] , 

( 86 ) 

(|pX-9« + <V) + o[3]. 

(87) 


The nonlinear dynamical equations are completed by the integrated Boltzmann multipole equations given in Section 
V - see Eq. dzl)- For neutrinos (£ > 2): 

0 = fi<^^> + |0n^^ + ^ 

- + (^ + 


■ N 

U 


(£-!)(£+!)(£ +2) 

(2£ + 3)(2£ + 5) ^ (2£ + 3) 

For photons, the quadrupole evolution equation is 




( 88 ) 


+ 40^a. ^ ^ 


= —Tlr^cr^ 


10 R 


' 0[3]. 


(89) 
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In the free-streaming case = 0, equation reduces to the result first given in |Q. This quadrupole evolution 
equation is central to the proof that almost-isotropy of the CMB after last scattering implies almost-homogeneity of 
the universe [^ . 

The higher multipoles {£ > 3) evolve according to 


\2it H" oj 

- + {£ + 
y2l. “h oj 


be 


(£-!)(£+!)(£ + 2 ) 

(2£ + 3)(2£+5) ' (2^ + 3) 




- Tip: (Trj-. 




1 


2£-f 3 


n 


Aia„ 


0[3]. 


(90) 


For £ = 3, the second term in square brackets on the right of Eq. must be multiplied by The temperature 
fluctuation multipoles are determined in principle from the radiation dynamic multipoles 11^^ via equations ( |7^ ) 
and d?^). 

These equations show in a transparent and explicitly 1+3 covariant and gauge-invariant form precisely which physical 
effects are directly responsible for the evolution of CMB anisotropies in an inhomogeneous universe. They show how 
the matter content of the universe generates anisotropies. This happens directly through direct interaction of matter 
with the radiation, as encoded in the Thomson scattering terms on the right of equations (|^, (1^), (|^ ) and (|^). 
And it happens indirectly, as matter generates inhomogeneities in the gravitational field via the field equations 
(0) and the evolution equation (Q) for the baryonic velocity u“. This in turn feeds back into the multipole equations 
via the kinematic quantities, the baryonic velocity z;“, and the spatial gradient in the dipole equation (|^). 

The coupling of the multipole equations themselves provides an up and down cascade of effects, shown in general by 
equation (p0[). Power is transmitted to the £-multipole by lower multipoles through the dominant (linear) distortion 
term ^ as well as through nonlinear terms coupled to the 4-acceleration baryonic velocity 

and shear Simultaneously, power cascades down from higher multipoles through the 

linear divergence term (divll)"^'^, and the nonlinear terms coupled to A“, and (Note that the vorticity coupling 
does not transmit across multipole levels.) 

The equations for the radiation (and neutrino) multipoles generalize the equations given by Challinor and Lasenby 
to which they reduce when we remove all terms 0{ev^) and O(e^). In this case, i.e. FLRW-linearization, there 
is major simplification of the equations: 


and for £ > 3 


Pn + +div%j 

C + + (div^J“ 

+ diJTrf + + Ig (divn)“' 


0 , 




(^ + 1 ) 

(2£-h3) 


(divn)^*^ 




(91) 

(92) 

(93) 


(94) 


These linearized equations, together with the linearized equations governing the kinematic and free gravitational 
quantities, given by equations (ra~(l^ with zero right hand sides, may be covariantly split into scalar, vector and 
tensor modes, as described in [|4^ , p5 ]^ ] . The modes can then be expanded in covariant eigentensors of the comoving 
Laplacian, and the Fourier coefficients obey ordinary differential equations, facilitating numerical integration. Such 
integrations are performed for scalar modes by Challinor and Lasenby p^, with further analytical results given in 

US 

However, in the nonlinear case, it is no longer possible to split into scalar, vector and tensor modes A 

simple illustration of this arises in dust spacetimes, which may be considered as a simplified model after last scattering 
if we neglect the dynamical effects of baryons, radiation and neutrinos. If one attempts to carry over the linearized 
scalar-mode conditions [|3|j36) 



— 0 — Tl^b 5 
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into the nonlinear regime, it turns out that a non-terminating chain of integrability conditions must be satisfied, so 
that the models are in general inconsistent unless they have high symmetry ||^,^. Thus, even in this simple case, it 
is not possible to isolate scalar modes. In particular, gravitational radiation, with cml Hab yf 0 (see must in 

general be present. 

The generalized equations given above can form the basis for investigating the implications of nonlinear dynamical 
effects in general, and second-order effects against an FLRW background in particular. More quantitative and detailed 
investigations along these lines are taken up in further papers. Here we will confine ourselves to a qualitative analysis. 


A. Nonlinear effects on kinematic, gravitational and dynamic quantities 


-ab 


Evolution of the expansion of the universe 0, given by equation ( 0 ) , is retarded by the nonlinear shear term —crabcr° 
and accelerated by the nonlinear vector terms +AaA°‘ and +2L0aUJ°' (see also Q). The vorticity evolution equation (||) 
has a nonlinear coupling aab^^ of vorticity to shear, whose effect will depend on the alignment of vorticity relative to 
the shear eigendirections. The shear evolution equation & has tensor-tensor and vector-vector type couplings, which 
are the tensor counterpart of similar terms in the expansion evolution. But in addition, relative velocity effects enter 
via the total anisotropic stress term. From equation (|8^), we see that baryonic and cold dark matter contributions 
of the form pv^°‘v^'> to the shear evolution arise at the nonlinear level. The constraint equations (|^) and (|^) show 
that acceleration and vorticity provide scalar (A“a;o) and vector ([w, nonlinear source terms for respectively the 
vorticity and shear. 

The free gravitational fields, which 1-1-3 covariantly describe tidal forces and gravitational radiation (see P,0,p5|-^), 
and therefore in particular control the tensor contribution to CMB anisotropies, are governed by the Maxwell-like 
equations dO). ® and (0)- This is the foundation for the electromagnetic analogy. The role of nonlinear 

coupling terms in these equations is more complicated - see for a full discussion. Here we note that nonlinear 
couplings of the shear and vorticity to the energy flux and gravito-magnetic field act as source terms for the gravito- 
electric field - see Eq. dil), while nonlinear couplings of the shear and vorticity to the anisotropic stress and 
gravito-electric field act as source terms for the gravito-magnetic field - see Eq. ©• 

From equations (^)-(Q), we see that for baryonic and cold dark matter, nonlinear relative velocity terms act as a 
source for the linear parts of the evolution equations for energy density and relative velocity. While the 4-acceleration 
Aa is involved in correction terms in all these equations, the vorticity and shear aab only enter nonlinear corrections 
of the velocity equations, and not the energy density equations. This reflects the fact that vorticity and shear are 
volume-preserving. The kinematic corrections to the evolution of matter relative velocity are of the form AaV°', [ta, v]a 
and cTabV^. For the massless species, as shown by equations dil), and (^^, (1^), the same form of corrections 


arises in the energy flux evolution, since energy flux is of the form when the photon and neutrino frames are 

chosen as the energy frame. Vorticity also does not affect energy density, but shear does, owing to the intrinsic 
anisotropic stress of photons and neutrinos, which couples with the shear. 

Baryonic and radiation conservation equations are both affected by nonlinear Thomson correction terms, which in¬ 
volve a coupling of the baryonic relative velocitv to the radiation energy density, momentum density and anisotropic 
stress. In particular, we note that there is a nonzero energy density transfer due to Thomson scattering at second-order. 


B. Nonlinear effects on radiation multipoles 

Nonlinear Thomson scattering corrections also affect the evolution of the radiation quadrupole 7r“^, as shown by 
equation (pl). In this case, the baryonic relative velocity couples to the radiation dipole and octopole Note 

also the correction to the linear Thomson term n^oyTr^^, in agreement with [ p9p^ ]. This correction arises from 
incorporating anisotropic effects in the scattering integral (while neglecting polarization effects, as noted earlier). 

The general evolution equation ( po|) for the radiation dynamic multipoles H^^ shows that five successive multipoles, 
i.e. for i — 2, ■■■,£+ 2, are linked together in the nonlinear case. Furthermore, the 4-acceleration Aa couples to the 
£±1 multipoles, the vorticity Wq couples to the £ multipole, and the shear Gab couples to the £±2 and £ multipoles. 
All of these couplings are nonlinear, except for ^ = 1 in the case of Ag, and ^ = 2 in the case of Gab- These latter 
couplings that survive linearization are shown in the dipole equation ( p7|) (i.e. PjiA°‘) and the quadrupole equation 
( p^ ) (i.e. /ORcr“^). The latter term drives Silk damping during the decoupling process Q. Nonlinear corrections 
introduce additional acceleration and shear terms. Vorticity corrections are purely nonlinear, i.e. vorticity has no 
direct effect at the linear level, and a linear approach could produce the false impression that vorticity has no direct 
effect at all on the evolution of CMB anisotropies. However, for very high £, i.e. on very small angular scales, the 
nonlinear vorticity term could in principle be non-negligible. 
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The disappearance of most of the kinematic terms upon linearization is further reflected in the fact that the 
linearized equations link only three successive moments, i.e. £, £ ± 1. This is clearly seen in equation (p4|). 

In addition to Aa and uja, there is a further vector coupling at the nonlinear level, i.e. the coupling of the baryonic 
velocity to the ^ ± 1 multipoles in the Thomson scattering source term of the evolution equation ( po|) . In fact these 
nonlinear velocity corrections are of precisely the same tensorial form as the acceleration corrections on the left hand 
side, only with different weighting factors. Linearization, by removing these terms, also has the effect of removing the 
nonlinear contribution of the radiation multipoles to the collision multipole 

One notable feature of the nonlinear terms is that some of them scale like £ for large £, as already noted in the 
case of vorticity. There are no purely linear terms with this property, which has an important consequence, i.e. that 
for very high £ multipoles (corresponding to very small angular scales in CMB observations), certain nonlinear terms 
can reach the same order of magnitude as the linear contributions. (Note that the same effect applies to the neutrino 
background.) The relevant nonlinear terms in Eq. ( ^0|) are (for £ ^ 1): 


The observable imprint of this effect will be made after last scattering. In the free-streaming era, it is reasonable to 
neglect the vorticity relative to the shear. We can remove the acceleration term by choosing as the dynamically 
dominant cold dark matter frame (i.e. choosing ?;“ = 0), as in It follows from equations (^) and (Q) that the 
nonlinear correction to the rate of change of the linearized temperature fluctuation multipoles is 




4.-4/ 


'(1 


^bc'^ 


bcA^ 




for ^ > 1. 


(95) 


The linear solutions for and aab can be used on the right hand side to estimate the correction to second order. Its 
effect on observed anisotropies will be estimated by integrating from last scattering to now. (See for the 

relation between the ta^ and the angular correlations). In the case of scalar perturbations, these solutions are given 
by Challinor and Lasenby |3^ (see also f^ ). 

Finally, we note that the well-known Vishniac and Rees-Sciama second-order effects also become significant at high 
£, and can eventually dominate the linear contributions to CMB anisotropies on small enough angular scales (typically 
£ > 10^ or more) ||. 


C. Temperature fluctuation multipoles 

We can normalize the radiation dynamic multipoles to define the dimensionless multipoles {£ > 1) 

Thus the are equal to the temperature fluctuation multipoles plus nonlinear corrections. In terms of these 
quantities, the hierarchy of radiation multipoles becomes: 


T 

T 


rj~a 


/j~ab 


-i0-iD,r“ 

- « - T“) + 0[3], 

-4 f ^ + 10 ) r“ - ^ - I« - T“) 


2^ .... 2 ^ ^ab , , ^la 8 rrdb 


+ - UbT'^>> + [a;,r]“ - |T“ 


T 


+ 0[3] 


_4 ( ^ + 10 ) - 0-“^ - 


ab 


T 




D,T 


-b -b C>[3] 


(96) 


(97) 


(98) 


and, for £ > 3: 
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'^Bb 


{£+!){£-2) ^ 
(2£ + 3) 


^ 2 £ + 3 

-{£ + 3)^(“^r^^-i> + + {£ + 2)cr<“^“^-iT^^-=> 

5£ 


(£-!){£+!){£ +2) 

(2£ + 3)(2£ + 5) (2e + 3) 




(99) 


For 1 = 3, the Thomson term must be multiplied by |. 

The nonlinear multipole equations given in this form show more clearly the evolution of temperature anisotropies 
(including the monopole, i.e. the average temperature T). Althought the only determine the actual temperature 
fluctuations to linear order, they are a useful dimensionless measure of anisotropy. Furthermore, equations (|^)- 
( p^ apply as the evolution equations for temperature fluctuation multipoles when the radiation anisotropy is small 
(i.e. = A‘), but the spacetime inhomogeneity and anisotropy are not restricted. This includes the particular 

case of small CMB anisotropies in general Bianchi universes, or in perturbed Bianchi universes. 

FLRW-linearization, i.e. the case when only first order effects relative to the FLRW limit are considered, reduces 
the above equations to: 


and, for £ > 3: 


T ^ 
T " 

i -40 - , 

( 100 ) 


, - 

( 101 ) 

fab p 

. -a^b _ J 3 {a^h> _ _ ^n^cr^T‘^b ^ 

( 102 ) 



(103) 


These are the 1+3 covariant and gauge-invariant multipole generalizations of the Fourier mode formulation of the 
i nteg ra ted Boltzmann equations used in the standard literature (see e.g. 0 and the references therein). Equations 
( |100D -(102) were given in @ in the free-streaming case = 0. 

As noted before, there is still a gauge freedom here associated with the choice of 4-velocity Given any physical 
choice for this 4-velocity which tends to the preferred 4-velocity in the FLRW limit, the .^ > 1 equations are gauge- 
invariant. 


VII. CONCLUSIONS 

We have used a covariant Lagrangian approach, in which all the relevant physical and geometric quantities occur 
directly and transparently, as PSTF tensors measured in the comoving rest space. There is no restriction on the 
deviation of geometric and physical quantities from FLRW limiting values, so that arbitrary nonlinear behavior may 
in principle be treated. We have derived the corresponding equations governing the generation and evolution of 
inhomogeneities and CMB anisotropies in nonlinear generality, without a priori restrictions on spacetime geometry or 
specific assumptions about early-universe particle physics, structure formation history, etc. Thus we have developed 
a useful approach to the analysis of local nonlinear effects in CMB anisotropies, with the clarity and transparency 
arising from 3+ 1 covariance. The equations are readily linearized in a gauge-invariant way, and then the methods of 
may be used to expand in scalar modes and regain well-known first-order results | 6 |] (see also f^ , ^ , ^ ). 

This approach allowed us to identify and qualitatively describe some of the key local nonlinear effects, and more 
quantitative results will be considered in further papers. We calculated the nonlinear form of Thomson scattering 
multipoles (given the initial simplifying assumption of no polarization), revealing the new effect of coupling between 
the baryonic bulk velocity and radiation brightness multipoles of order £ ± 1. We also found the nonlinear effects of 
relative velocities of particle species on the dynamic quantities that source the gravitational field. These effects also 
operate on the conservation equations, including evolution equations for the relative velocities of baryonic and cold 
dark matter. 

Nonlinear effects come together in the hierarchy of evolution equations for the radiation dynamic (brightness) 
multipoles, which determine the CMB temperature anisotropies. In addition to the nonlinear Thomson contribution, 
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we identified nonlinear couplings of the kinematic quantities to the multipoles of order ^ ± 2, £ ± 1, and £. These 
quantities themselves are governed by nonlinear evolution equations, which provides part of the link between CMB 
anisotropies and inhomogeneities in the gravitational field and sources. The link is also carried by the spatial gradient 
of radiation energy density (equivalently, average radiation all-sky temperature), and the baryonic relative velocity. 
Furthermore, there is internal up- and down-transmission of power within the multipole hierarchy, supported by the 
kinematic couplings as well as by distortion and divergence derivatives of the multipoles. 

We used our analysis of the radiation multipoles to identify new effects that operate at high £. In particular, we 
showed that there is a nonlinear shear correction effect on small angular scales, whose impact on the angular power 
spectrum was qualitatively described. The quantitative analysis of this and other nonlinear effects is a subject of 
further research. 
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APPENDIX A: EXACT NONLINEAR RELATIVE VELOCITY EQUATIONS 

Change in 4-velocity: 

Ua=l{ua + Va) where 7 = (I, WqU® = 0. (Al) 

Change in fundamental algebraic tensors: 

Kb = Kb + 7 ^ [v'^UaUb + ‘2U(^aVb) + VaVb] , (A 2 ) 

^abc 'y^abc T T T • (-^^) 

Transformed kinematic quantities are defined by 

^ b'^a — 3 T ^ab T ^abc^ -^a^b ; 

which implies, using Va 7 = j^v^VaVb, and Eq. (^), the following kinematic transformations [^ : 

0 = 7 © - 1 - 7 (div V + AKa) + , (A4) 

Aa = 'J^Aa + 7^ {U(a) + ^QVa + (JabV^ - [w, v]a + + A^Vb + abcV^v‘') Ua 

+ ^{divv)Va + i['U,CUrlu]a -b U*'D(bUa)} -b j‘^W{Ua + Va) , (A5) 

Wo = 7^ { (l - Wo - ^curl?;a -b \vb (2w'' - curl?;^) Ua -b \vbLO^Va 

-b \ [A,v\a -b \ [v,v\a + \eabc(^''dV^'v’^) , (A 6 ) 

dab = ICFab + 7(1 + l'^)u{a<yb)cK -b 7^^(a [vb) + v'^Ub)] 

+ l^{aVb) - \Kb [AcV^ -b 7 ^ (W - VcV’')] 

+ j^UaUb -b - v^t)‘*D(cUd) -b ( 7 ^ - “ l) W] 

+ l^U^aVb) [AcV'' + (TcdV^'v'^ - Uc'c" -b ( 7 ^ - i) W] 

+ \l'^VaVb [div?; - AaV^' + 7 ^ ( 87 ^ - l) W] -b J^Vi^aVb) + V^l'^'U(aV(b)) 

+ 7^^^(aO-b)c'C° - 7 ^[w, w](o {ub) -b v'^Ub)) + 2'-/^ v'^D (^a.V (a) {l^b) + Mb) } , (A7) 

where 

W = VcV^ -b |?;^divr; -b u'^u'^D^cUd) . 

Transformed dynamic quantities are [^2|: 
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P = P + 7^ K(p +p) - 2 gaf“ + TTabU^U^] , 

(AS) 

p = p + ^ 7 ^ [u^(p + p) - 2qaV‘^ + 7ra6u“u^] , 

Qa = 79o - JT^abV^ “ 7^ [(P + P) “ + TTbcV^V^"] Va 

(A9) 

- 7 ^ [u^(p+p) - (1 + v^)qbv'' + TTbcV^v'^] Ua , 

^ab = T^ab + 27^U%e(a [ub) + Q{aUb) “ ‘^I’^Qla^b) 

[v^{p + p) + hab 

+ ^ 7 "^ [2v^{p + p) - Av'^qaV‘^ + (3 - V^)TTadV''v'^] UaUb 
+ § 7 "^ [2v^{p + p) - (1 + 3v‘^)qcV‘' + 2TradV‘'v‘^] U(^aVb) 

(AlO) 


+ [(3 - v‘^){p + p) - + 2TTcdV''v'^] VaVb . (All) 

gravito-electric/magnetic field: using 

= 4 [u^aU^^ + + 2eabeU^^H<^^^ + 

= 4 Eb]‘'^ + 2eabeU^^H‘‘^‘^ + , 

we find the transformation ||^ ]: 

Eab = 7^ {(1 + V^)Eab + \p‘£cd{aHb)'^ + 2Ac(aMb) 

+ (UaMb + hab)EcdV‘^ - 2Ac(aUb) + 2u^aSb)cdH‘^''Ve] } , (A12) 

Hab = 7^ {(1 + V^)Hab + \—‘^£cd(aEb)'^ + ‘^Hc(aUb) 

+ {UaUb + hab)HcdV‘^ - 2iJc(a^^b) “ ‘2U{aSb)cdE‘^‘'Ve] } ■ (A13) 

This may be compared with the electromagnetic transformation 

Ea=l {Ea + [v,H]a + v'^EbUa} , 

Ha=l {Ha - [U, E]a + v'" HbUa) , 

where 

Eab — ‘^'^[aEb] “t“ ^ahcH — 2'll]^aEb'\ T ^abcH . 

Note that all the transformations above are given explicitly in terms of irreducible quantities (i.e. irreducible in the 
original M“-frame). 
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